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Abstract 

We analyse the effects of the propagating Alfven waves, arising due to non-zero mean magnetic 
fields, on the nonequilibrium steady states of three-dimensional (3d) homogeneous Magnetohydro- 
dynamic (MHD) turbulence. In particular, the effects of Alfven waves on the universal properties of 
3dMHD turbulence are studied in a one-loop self-consistent mode-coupling approach. We calculate 
the kinetic- and magnetic energy-spectra. We find that even in the presence of a mean magnetic 
field the energy spectra are Kolmogorov-like, i.e., scale as /c _5//3 in the inertial range where k is a 
Fourier wavevector belonging to the inertial range. We also elucidate the multiscaling of the struc- 
ture functions in a log-normal model by evaluating the relevant intermittency exponents, and our 
results suggest that the multiscaling deviations from the simple Kolmogorov scaling of the struc- 
ture functions decrease with increasing strength of the mean magnetic field. Our results compare 
favourably with many existing numerical and observational results. 

PACS numbers: 64.60Ht,0.5.70.Ln,47.65+a 



I. INTRODUCTION 



The effects of propagating waves on the statistical properties of systems out of equilib- 
rium remain an important topic of discussion. In the context of a coupled spin model in 
one-dimension (d = 1) l| it has been shown that the presence of such waves leads to weak 
dynamic scaling in that model. In contrast, in a coupled Burgers-like model in d = 1 prop- 
agating waves do not affect the scaling properties of the correlation functions at all [2]. So 
far such issues have been considered only within very simplified one-dimensional nonequilib- 
rium models . Magnetohydrodynamic (3dMHD) turbulence, which is a hydrodynamic 
description of the coupled evolution of the velocity fields u and the magnetic fields b in a 
quasi-neutral plasma, stands as a very good candidate for a natural system with propagating 
waves in three dimensions as most of its natural realizations have proagating Alfven waves 
arising due to the presence of mean magnetic fields. Examples of such physical situations 
include solar wind, neutral plasma in fusion confinement devices etc. The presence of the 
propagating Alfven modes, in addition to the usual dissipative modes due to the fluid and 
magnetic viscosities makes MHD turbulence a good natural example to study the interplay 
between the propagating and the dissipative modes in a system and their combined effects 
on the scaling properties of the correlation and structure functions, which are important 
issues from the point of view of nonequilibrium statistical mechanics. 

The scaling of magnetic- [-£/(,(&)] and kinetic- [E u (k)} energy spectra in the inertial range 
(i.e., wavevector k lies in the region L^ 1 <C k <C 1]^, L and rjo being the integral scale 
given by the system size and the dissipation scale, respectively) in 3<iMHD in the presence 
of Alfven waves originating due to a mean magnetic field B Q remains controversial till the 
date. Numerical simulations, due to the lack of sufficient resolutions failed to conclusively 
distinguish between the Kolmogorov's and Kraichnan's predictions (see below). In this 
communication, within one-loop self-consistent mode-coupling (SCMC) approximations, we 
obtain the following results for homogeneous but anisotropic (due to the mean magnetic 
field) 3dMHD: 

• The bare Alfven wave speed, proportional to B Q , renormalizes to acquire a singular 
/c-dependence to become ~ B Q k~ 1 ^ where k is a Fourier wavevector belonging to the 
inertial (scaling) range. From this result we are able to conclude that even in the 
presence of a mean magnetic field the kinetic- and the magnetic- energy spectra scale 
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as k 5//3 in the inertial range, identical with the situation without any mean magnetic 
field. 



• The dimensionless Kolmogorov's constants for the kinetic- and magnetic-energy spec- 
tra depends on a dimensionless parameter (3 which we identify as the ratio of the 
renormalised Alfven wave speed and the renormalised viscosities (see below). 

• The intermittency exponents of the Elsasser fields = u ± b which approximately 
and qualitatively characterise the multiscaling properties of the structure functions in 
a log-normal model depend on the parameter j3 and decrease with increasing /3. 

Thus our results show that Alfven waves in 3<iMHD do not affect the scaling properties 
of the two-point correlation functions. However, the multiscaling properties of the structure 
functions are shown to be affected by the mean magnetic fields. 

The famous Kolmogorov's arguments j^J for fluid turbulence can be easily extended to 
3<iMHD turbulence. In MHD, in the unit where mass density p = 1, the kinetic energy 
dissipation rate per unit mass (ex) and the magnetic energy dissipation rate per unit mass 
(ejvf) have same physical dimensions. Thus, as in fluid turbulence, by claiming that the 
structure functions of the velocity and magnetic field differences in the inertial range must 
be constructed out of the mean energy dissipation rate e per unit mass and the local length 
scale r (belonging to the inertial range) one obtains for the n-th order structure function 

DO 

S«(r) = ([oi(x + r) - a^x)]*]") ~ (er)"/ 3 , rj D « r « L (1) 

where e = ekOtem and a = u,b for the velocity and the magnetic fields, tjd is the (small) 
dissipation scale and L is the (large) system size. This yields, for the energy spectra in the 
inertial range (as a function of wavevector k) 

E a {k) ~ k- 5/3 , a = u,b. (2) 

This is known as the K41 theory in the relevant literature. However, in the presence of a 
mean magnetic field B Q , the Alfven waves are generated (see below) with the propagation 



speed ~ B a . Thus in such a case in addition to the usual Kolmogorov time scale ~ k 2 ^ 3 
there exists another time scale constructed out of the mean magnetic field B Q , known as the 
Kraichnan time scale ~ (Bok) -1 Kraichnan argued that this time scale would determine 
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the energy cascade process and hence would enter in the expression of the structure function 
yielding E a {k) ~ k~ 3 ^ 2 , a — u,b. There has not been any satisfactory resolution of this issue 
till this date; due to the particularly difficult vectorial nature of the 3<iMHD equations (see 
below) it is rather difficult to achieve high Reynolds number in Direct Numerical Solutions 
(DNS) of the 3<iMHD equations. Numerical solutions of an MHD shell-model in the presence 
of a small mean magnetic field-like term did not find any dependence of the multiscaling of 
the structure functions on the mean magnetic field p . Analytically, it has been shown within 
the context of a Id coupled Burgers-like model that the energy spectra are independent 
of a mean magnetic field and in case of Id MHD turbulence it exhibits the Kolmogorov scal- 
ing for the energy spectra. Similar conclusions followed from an analogous one-dimensional 
model 0. Various phenomenological approaches, including weak turbulence theories and 
three-wave interaction models, yield, in general, mean magnetic field dependences of the 

when the mean magnetic field is strong. However, these the- 



energy spectra 
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ories, despite their predictions are either not directly derivable from the underlying 3<iMHD 
Eqs. of motion or involve additional assumptions on the flow fields. Analyses starting from 
the 3<iMHD Eqs. in this regard are still lacking. Most observational results on astrophysical 
systems seem to favour the K41 results jl^ . Simulations of incompressible MHD 15| find 
results close to the K41 result. Simulations of compressible MHD |l6j], even though in some 
cases find energy spectra closer to K41, in general yields a less clear picture. Recent numer- 

in presence of finite magnetic helicity structure 
functions parallel and perpendicular to the mean magnetic fields are affected differently by 
the mean magnetic field. In this paper, we address some of these issues by starting from 
the 3c/MHD Eqs. without making any further assumptions on the velocity and the magnetic 
fields, except for the validity of the perturbative approaches. In particular, we show by 
applying one-loop mode coupling methods on the 3c/MHD equations with a mean magnetic 
field and in the absence of any magnetic helicity that the one-dimensional energy spectra 
in 3(iMHD turbulence are independent of the mean magnetic field B Q and scale as k~ 5 ^ 3 in 
the inertial range where k is a Fourier wavevector belonging to the inertial range. We then 
proceed to calculate the Kolmogorov's constants for the kinetic- and the magnetic-energy 
spectra and show that they depend on B Q . Lastly, we calculate the intermittency exponents 
for the velocity and the magnetic fields and find that they decrease with increasing B . We 
do not distinguish between the longitudinal and the transverse structure functions. Even 
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though for analytical convenience we assume a weak mean magnetic field, we are able to 
obtain new and interesting results concerning scaling and multiscaling in 3dMHD in the 
presence of a mean magnetic field of small magnitude. In this respect our results can be 
considered as complementary to some of the existing results The rest of the 

paper is organised as follows: In SecjHJwe discuss the 3g?MHD equations for incompressible 
fluids. In Sec]V]we show that for incompressible fluids the bare Alfven wave speed ~ B Q 
renormalizes to pick a correction fc _1//3 in the inertial range. This, as we argue, implies 
that the energy spectra scale as k~ 5 ^ 3 in the inertial range even in the presence of a mean 
magnetic field. In Sec I VII we calculate the Kolmogorov's constants for the kinetic and the 
magnetic energy spectra. We introduce a parameter (3 which is the dimensionless ratio of 
the renormalised mean magnetic field and renormalised viscosities and show that the Kol- 
mogorov's constants depend on (3. In Sec JVIII we elucidate the multiscaling properties of the 
structure functions by calculating the intermittency exponents and show that they decrease 
with increasing (3, i.e., with increasing B a . This suggests that a mean magnetic field tends 
to reduce multiscaling corrections to the K41 results for the structure functions. In Sec JVIIll 
we summarize our results. 



II. MODEL EQUATIONS 

We begin by writing down the 3<iMHD equations for the velocity fields u and the magnetic 
fields b: The velocity field u is governed by the Navier-Stokes equation modified by the 
inclusion of the Lorentz force ^| 

du / v~~r \ Vp (V x b) x b „ 

— + u-Vu = ^ + ^ '- + z/V 2 u + f, 3 

ot p p 



and the dynamics of the magnetic field b is governed by the Induction equation [18( con- 
structed out of the Ohm's law for a moving frame and the Ampere's law: 

db 

— = V x (u x b) + r/V 2 b + g. (4) 

Here, p is the mass density, p is the pressure, v is the fluid viscosity and rj is the magnetic vis- 
cosity (inversely proportional to the electrical conductivity of the fluid medium concerned). 
Functions f and g are external forces needed to maintain a statistical steady state. In the 
present approach these are taken to be stochastic forces. We assume them to be zero-mean 
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and Gaussian distributed with specified variances (see below). In addition to Eqs. Q and 
(jlj) we also have V • b = (Maxwell's equation) and, for incompressible fluids V ■ u = 0. 

If the magnetic fields b(x, t) are such that (b) = B Q (a constant vector) then replacing 
b by b + B Q where now (b) = 0, in Eqs. flHJ) and (@J) one obtains additional linear terms 
proportional to wavevector k leading to wave-like excitations, known as Alfven waves. The 
resulting equations are 

du x . ^ . Vp , (Vxb)xb (V x b) x B „ 

— +Ai(u-Vu = - + X 2 - + - - + z/V 2 u + f, 5 

ot p p p 

and 

db 

— = A 3 V x (u x b) + V x (u x B ) + r/V 2 b + g. (6) 

with (b(x, t)) = 0. The parameters Ai, A2, A3 are kept for book keeping purposes and can 
be set to unity (see SeclIV|). Note that on dropping the nonlinear and the dissipative terms 
from Eqs. © and (jOJ) the resulting linear coupled partial differential equations admit wave- 
like solutions with dispersion relation linear in wavevector k. These are known as the Alfven 
waves ^| in the literature which propagate with speed proportional to B . In the following 
sections we would calculate the kinetic and the magnetic energy spectra, the Kolmogorov's 
constants and the intermittency exponents in the presence of the Alfven waves, i.e., for 



III. CORRELATION AND STRUCTURE FUNCTIONS IN MHD 

In the statistical steady state the time dependent correlation functions of u and b exhibit 
scaling which are characterized by the roughness exponents Xu and Xb respectively, and the 
dynamic exponent z. In terms of the scaling exponents Xu, Xb and z the velocity and the 
magnetic field correlators have the form (as a function of wavevector k and frequency to) 

C%{k,u) = ( Ui {k,u) Uj {-k, -u)) = DoPyk-^-'Mk'/u), 
C^(k,u) = (&,(k,^(-k, -a;)) = D b P ij k- d - 2 x>- z f b (k t /u), (7) 

where f u and /& are scaling functions, Py is the transverse projection operator: Py = 
8{j — kikj/k 2 to account for the incompressibility of the fields. The kinetic- and magnetic- 
energy are simply related to the correlation functions: They are just the equal time velocity 
and magnetic field auto correlation functions multiplied by appropriate phasefactors. We 
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now set out to find whether the Alfven waves are relevant perturbations on the system. It 
should be noted that the correlators (J7J) are chosen as they would be in the fully isotropic case. 
In the presence of a mean magnetic field B Q there would however be additional anisotropic 
terms in the expression for the correlators above. In the small B Q limit the lowest order 
corrections to the expressions (JJJ) are 0(B ) 2 . Neglecting these correction terms does not 
affect our scaling analyses below, and in any case, we are interested to see whether anisotropic 
Alfven waves are relevant perturbations in the large scale, long time limit on the isotropic 
3<iMHD as characterised by expressions (J7J). Therefore, it suffices for us to work with the 
expressions (J2J) for our mode-coupling analyses below. 

A complete characterisation of the nonequilibrium steady state (NESS) of MHD, however, 
requires informations about, not just the scaling of the energy-spectra, but also of the n-th 
order equal time structure functions of the velocity and the magnetic fields in the inertial 
range. These are defined by 

SZ(r) = (\[ai(x + T)-a4(x)]ri\ n ), a = u,b, r] D < r < L, (8) 

which scale as r^™ where r belongs to the inertial range (rjjj C r < L). According to the 
Kolmogorov theory (K41) the multiscaling exponents £™ = n /3; i>e., they are linear in n. 
Subsequent numerical and observational results for homogeneous and isotropic 3dMHD, i.e., 



without any mean magnetic fields, (2Q, |2l| su gge sted deviations from the K41 results which 



are similar to those of pure fluid turbulence 



22[. Much less results for the multiscaling of 



the structure functions are available when there are mean magnetic fields. Recent numerical 
simulations of Miiller et al suggested that structure functions parallel and perpendicu- 
lar to the mean magnetic fields are differently affected by it when there is a finite magnetic 
helicity. Below we investigate the issue of multiscaling in presence of a mean magnetic field 
B within the context of a log-normal model by evaluating the relevant intermittency expo- 
nents which are found to depend explicitly on B . In our analyses we ignore the distinction 
between the structure functions parallel and perpendicular to the direction of the mean 
magnetic field; nevertheless, as we discuss below, our results are significantly new. 
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IV. SYMMETRIES OF THE EQUATIONS OF MOTION 



We begin by re-expressing 3<iMHD Eqs. © and (HJ) by writing the magnetic fields b as 
a sum of a space-time dependent part and a constant vector: b(x, t) —>■ b(x, t) + B Q . In 
terms of these fields and parameters, the 3<iMHD Eqs. become 

du _ . . Vp , (Vxb)xb (V x b) x B _„ „ 

— + Ai(u-V)u) = - + A 2 - + - - + zA7 2 u + f, (9) 

ot p p p 



and 



dh 

— = A 3 V x (u x b) + V x (u x B ) + r]V 2 b + g. (10) 



It should be noted that in Eqs. (JHI) and B Q is not the mean magnetic field in general; 
it would be so only if (b(x, £)) is zero in Eqs. Q and (jlOj) . Note that our above way of 
splitting the magnetic fields fields does not change the actual mean magnetic field in the 
system: It is still given by B Q + (b) = B Q . The Eqs. of motion Q and (fTU|) are invariant 
under the following continuous transformations |sl U^: 

re^Galilean transformation (TI): u(x, t) — > u(x + u Q t, t) + Uo, J^ — uo.V, and b — » h 
3l wi 




24, 



with Ai = A3 = 1 in Eqs. (0) and ©• This implies non-renormalization of 



• The transformation (Til) B — > B + A2<5, b(x, t) — > b(x, t) — 5, u — >• u. This allows 
one to work with the effective magnetic fields defined by V^b such that the coefficient 
of the Lorentz force vertex constructed out of the effective magnetic fields does not 
renormalize. This, therefore, ensures A2 can be set to unity P, |^ by treating all 
magnetic fields as effective fields. Here the shift 5 is a vector. A transformation 
similar to Til above exists in a problem of passive scalar turbulence j^ . 

The transformation Til essentially signifies the freedom to split the total magnetic fields 
as a sum of a constant part and a space time dependent part. It should be noted that the 
transformation Til keeps the mean magnetic field unchanged. In fact, nonrenormalization 
of A2 can be shown in a simpler way: Let us assume that under mode eliminations and 
rescaling A2 — > aA2- This scale factor of a can now be absorbed by redefining the units of 
the magnetic fields by y/ah — > h. Therefore, A 2 can be set to unity if all magnetic fields are 
considered as effective or rescaled magnetic fields. Since the Induction Eq. (JU is linear in 
the magnetic fields b(x, t), such rescaling leaves every conclusion unchanged. Of course, the 
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external force g is also scaled by a factor y/a which does not affect out analysis here as the 
assignment of canonical dimensions to various fields and parameters is done after absorbing 
A2 in the definition of b. More specifically under the rescaling x — > Zx, t — > l z t, U{ — > 
l Xu Ui, hi — > l Xb bi the bare parameters scale as Ai )3 (Z) — > / Xu+2:_1 Ai i 3, X 2 (l) — > l 2xb ~ Xv + z ~ 1 \2 
and _B — > / Xf, ~ Xll+z_1 I? . Since there are no fluctuation corrections to the nonlinearities 
-^l) -^2) A3, which are the consequences of the invariances under the transformations TI and 
Til, they can be kept invariant under rescaling of space and time as mentioned above leading 
to Xu = Xt = X an d X + z — 1- Thus, under naive rescaling the bare parameter B a ~ B Q l z ~ l . 
Furthermore, the invariance under the transformation Til and the resulting Ward identity 
ensure that different ways of implementing one-loop RG by having different values for (b) 
and B Q in Eqs. (jSJ) and (|TU|) subject to the same bare mean magnetic field B Q = B Q + (b). 
Since any renormalisation scheme must respect a freedom of choice as represented by the 
transformation Til, respective terms must scale the same way under the rescaling of space 
and time j^. Hence, due to fluctuation corrections B (l) ~ l x . Furthermore, since both 
B and b have the same scaling dimensions given by the exponent Xi t ne physical mean 
magnetic field in the system B Q = B Q + (b(x, t)), if it receives fluctuation corrections under 
mode eliminations, must scale as B Q (l) ~ l x , or if there are no fluctuation corrections to 
it due to some special symmetries (as in the one- dimensional Burgers-like model for MHD 
in Ref. P]) will be irrelevant (in an RG) sense as it will flow to zero as l z ~ l (since z for 
fully developed turbulence is less than unity). This clearly suggests the possibility of a 
renormalisation of the Alfven wave speed, akin to the renormalisation of the sound speed 

n 

in compressible fluid turbulence [221]. As we will see below, for 3dMHD there are infra-red 
singular fluctuation corrections to the bare mean magnetic field leading it to renormalise in a 
way consistent with our predictions from the Ward identities above. It should be noted that 
our analyses above is independent of the strength of the bare value of the mean magnetic 
field. Although, in the discussion above we have rescaled space isotropically (i.e., x,y and 
z coordinates are rescaled the same way) and analyze the scale dependence of the resulting 
effective parameters and the fields, it should be noted that such a rescaling does not imply 
that the effective parameters have isotropic structures. 
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V. ENERGY SPECTRA FOR INCOMPRESSIBLE FLUIDS 



We begin by writing down the 3<iMHD equations in the incompressible limit (i.e., V • u = 
0). We write down the equations of motion in k-space in terms of the Elsasser variables 
z ± = u ± b. The equations are (we take the mean magnetic field B Q to be along the z 
direction) 

d^ t M) - iB k z zf + iPipks Z 7(<1, t)z}(k -q,t) + V °+k 2 z+ + if_k 2 z^ = 9+(k, t), 
dZl t] + iB Q k z z- + iP lv k s J2 4(q. t)z~ (k - q, t) + V°+k 2 zi + rj°_k 2 z+ = 0f(k,t)(ll) 

The stochastic forces 9 l are linear combinations of fi, gi, the tensor Pi p is the transverse 
projection operator which appears due to the divergence-free conditions on u and b, and 

T)± — V ± 77. 

In the absence of any crosscorrelations between the velocity and the magnetic fields, the 
simplest choice for the noise variances, consistent with the divergence-free conditions on the 
velocity and the magnetic fields, are 

(0±(k,t)0±(-k,t)) = 2P lm D x k- y 5{t), 

(6?{k,t)0*{-k,t)) = 2P lm D 2 k-y5(t). (12) 

In Eqs. fll2|) we choose y > 0. In particular, the choice of y = d in <i-space dimensions 
ensures that the energy flux in the inertial range is a constant (up to a weak logarithmic 
dependence on wavevector k) which forms the basis of the K41 theory j^], and as a result, the 
parameters D\ and D2 pick up dimensions of energy dissipation rate per unit mass. Hence, 
we will use y = d in our calculations below. In experimental realisations of MHD, external 
forces act on the large scales. In analytical approaches, such forces are replaced by stochastic 
forces with variances given by Eqs. (fT2"|) for calculational convenience. It should however be 
noted that in numerical solutions of the Navier-Stokes equation driven by stochastic forces 
with variances similar to Eqs. (|T2*|) structure functions of the velocity fields are shown to 
exhibit multiscaling similar to the experimental results [23]. Our preliminary results from the 
numerical solutions of the isotropic 3<iMHD equations (i.e., no mean magnetic fields) yield 
multiscaling similar to those obtained from the 3dMHD Eqs. driven by large-scale forces 
I29I ] . Although no such numerical results exist for the case with a mean magnetic field, it will 
presumably be true for such a situation also. In a stochastically driven Langevin description, 
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correlation functions are proportional to the noise correlations and hence it is necessary to 
force the Induction Eq. © stochastically, as is common in the literature (see, e.g., 3, Q 
and references therein). For such noises in the absence of any mean magnetic field (i.e., for 
the isotropic case) it has been shown by using a one-loop mode coupling theory that the 
scaling exponents have values x = 1/3, z = 2/3 jq. 

The Eqs. of motion (|11|) are coupled even at the linear level leading to a bare propagator 
matrix for the Eqs. (fTTj) of the form (a V refers to bare quantities) 



We use a one-loop mode coupling theory which is conveniently formulated in terms of 
the self-energy matrix £ and the correlation functions given by Eqs. (J7|). The self-energy 
matrix £ is defined by G^ 1 = G" 1 — X where G is the renormalised propagator matrix. In 
terms of the scaling exponents \ and z, and the renormalized parameters the self-energy 
matrix £ is given by 



where B(k) is the renormalised or the effective Alfven wave speed. If there are diagrammatic 
corrections to the imaginary parts of Y,(k,uj) at frequency u = which are singular in the 
infra-red limit, then B(k) is different from B Q , the bare Alfven wave speed, else they are the 
same. 

There are two one-loop diagrams which contribute to the fluctuation corrections to 
E u (fc,w = 0) = —iB(k)k z + i] + k z . These are shown in Fig.(|T|). These have both real 
and imaginary parts at frequency to = 0. Thus there are fluctuation corrections to the bare 
Alfven wave speed which are singular in the small wavenumber limit. We assume, in the 
spirit of mode-coupling methods, B{k) = Bk~ s (we assume an isotropic scale dependence for 
B (k) for simplicity which suffices our purpose of finding the scale dependence of the effective 
Alfven wave speed). Clearly, if — s + 1 < z, the small wavenumber limit of the problem is 
dominated by underdamped waves, if — s + 1 > z the Alfven waves are damped out in the 
small wavenumber limit. In contrast if — s + 1 = z then in the small wavenumber limit both 
the propagating and the dissipative modes are present. Our analyses in Sec II VI clearly sug- 
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gest that if there are fluctuation corrections to B Q then the effective scale dependent Alfven 
wave speed B{1) ~ k~ x yielding s — %. Furthermore, since % + z = 1 we have — s + 1 = z 
leading to the co-existence of the underdamped Alfven waves and the dissipative modes in 
the large-scale, long-time limit. Note that this situation allows us to define a dimensionless 
parameter (3 = B/r] + where B and r] + are the renormalised Alfven speed and the viscosity 
respectively. 

In the SCMC approach vertex corrections are neglected. Lack of vertex renormalisations 
in the zero wavevector limit in 3c?MHD allows SCMC to yield exact relations between the 
scaling exponents x an d z as in the noisy Burger s/Kardar-Parisi- Zhang equation [2^|. In 
this model the nonlinearities and the noise variances do not renormalise, thus leading to 
z = 2/3, x — 1/3, satisfying the exponent identity x + z — 1- in t ne present case, due to 
the singular nature (in the small wavevector limit) of the bare noise correlations (|12p. they 
do not pick up any further singular corrections to their scaling; however the amplitudes get 
modified (this is similar to the results in Ref.js]). We denote the renormalized amplitudes 
by D and D, respectively. Furthermore, for the 3<iMHD Eqs. (J5J) and (jOJ) there are infra-red 
singular fluctuation corrections to B a [see the one-loop diagrams in Fig. (P)] leading to 
B (l) ~ l x consistent with our arguments above. The SCMC approach involves consistency 
in the scaling and the amplitudes of the mode coupling equations. It should be noted 
that the exponent values z = 2/3, x = 1/3 satisfy the mode coupling integral equations 
regardless of the strength of the mean magnetic field. For amplitude consistency the one- 
loop integrals are required to be evaluated. Due their complicated structure, we evaluated 
them by assuming that the strength of the mean magnetic field is small. Therefore, only 
our amplitude relations and not the scaling exponents are affected by the approximation of 
small mean magnetic fields. 

Demanding amplitude consistency in the mode coupling equations we obtain, 



V+ 

-r 



D 
B 

7] + 



B = B, 



d d(d + 2) 
D 



+ 



D _ (3 2 D 
d{d+2)r]l ~ 



2 

1--4 



(3 2 D 



;i + r 2 )(i 



°d{d + 2)rf + 

d ~ d(d+2)> 



-0.5r(l 2) +jPl (/3) 



l + r (i + r 2 )(i 



2 _i 2 ■ 
d d(d+2) - 



+ 0.5r(l - 2) + F 2 (f3) ' 



where Fi(/3) = 2(3 2 (1 + 3T 2 )(-l/d + 



i 



d(d+2) d(d+2)(d+4). 



d d{d + 2) r] 2 + d{d + 2)(d + 4) 



(13) 



F 2 (/3) = 2(3 2 (3 + T 2 )(-l/d + 
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+ 

z 



z + / \ Z + 




(b) 

FIG. 1: One-loop diagrams contributing to renormalization of . A single line refers to a propa- 
gator and a line with a filled circle refers to a correlator. 

d(d+2) ~ d(d+2)(d+4) ) are ^ ne /^-dependent parts of the amplitude-ratio. While calculating 
the amplitude consistency relations we worked in the limit of small (3 = B/rj + and the 
renormalised magnetic Prandtl number (2 = ^+~^~ ) close to unity. Away from these 
limits the amplitudes of the underlying one-loop integrals become much more complicated 
functions of (3 and also of the renormalised magnetic Prandtl number, but the qualitative 
picture remains unchanged. The main physical picture that emerges from the expressions 
(fTHJ) is that in the absence of a bare mean magnetic field B Q the effective Alfven wave speed 
is also zero which is a restatement of the fact that if the original theory is isotropic, so 
will be the renormalised theory. Moreover the renormalised Alfven speed increases with 
increasing B . With these results we are now in a position to calculate the energy spectra 
in the inertial range. We use the form of the effective (scale dependent) viscosity and the 
Alfven wave speed to obtain the equal time correlation function of zf in the long wavelength 
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limit. We define (in 3d) 

2Dk~ 3 P ij (k) 



C± + (k,u) = (zt(k,u)zi(-k,-co)) 



(u - Bk-^k,) 2 + rf + k A / 3 ' 



Crr {k ,u) = (z-{k,u)zr{-k,-u))- ~ L " W 



C+-(k,U>) = (z+(k, u) Z T(-k, -U)) = H ^ _ Bk _ 1/3kz) + 2 ^ + p/3]^j + Bk -l/3 kz) + v+k i)tf 

As discussed before, we have omitted anisotropic corrections of 0(B o ) 2 or of 0(/3) 2 to the 
correlation functions as we are trying to find out the relevance (in an RG sense) of Alfven 
waves on the scaling of the isotropic correlation functions in the absence of any mean mag- 
netic fields. Therefore, the equal time correlation functions have the following form (in 
d = 3): 

C±+(k,t = 0) = Crr(k,t = 0) = Dk' 3 " 2 ' 3 ^). (15) 

Therefore, the equal time autocorrelation functions of are independent of any mean 
magnetic field. This holds true regardless of the scaling of the noise variances. The equal 
time cross correlation function C^~(k,t = 0) requires more careful considerations. On 
integrating C^~(k,u) over all frequency lo one obtains (in 3d) 

a form which is valid in the inertial range. It is clear from the expression (JTfij) that both the 
real and the imaginary parts of C^~(k, t = 0) scale as k~ 3 ~ 2 / 3 in the inertial range. Thus the 
one-dimensional kinetic- and the magnetic-energy spectra (which are simply related to the 
correlators defined in Eqs. (|14p) scale as k~ 5 ^ 3 in the inertial range. The emerging physical 
picture is as follows: We find from the expression for C^~(k, uj) above that this, as a function 
of frequency u, has maxima at uj = ±Bk~ 1 / 3 k z and the width at half-maxima ~ k 2 / 3 . In 
contrast, the auto correlations of z + ,z _ are maximum at uj = and their widths scale as 
~ k 2 l 3 . Thus in the long wavelength limit the width and the location of the maxima of 
C^~(k,u) scale in the same way leading to the presence of the underdamped Alfven waves 
in the hydrodynamic limit. Therefore, it immediately follows that the kinetic- and the 
magnetic-energy spectra, being linear combinations of the correlators discussed above times 
appropriate phase factors, scale as k~ 5 ^ 3 even in the presence of a non-zero mean magnetic 
field. It should be noted that we have used a small (3 approximation to arrive at our results for 
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the self-consistent amplitude relations. For a finite (3 one would require to work with a fully 
anisotropic form of the correlation functions and obtain self-consistent relations for scaling 
and anisotropic amplitudes, a task, which is analytically challenging, remains to be done in 
the future. However, based on our calculations above, the exponent identity = 1 an d the 
Ward identity suggesting that the renormalised Alfven wave speed should scale as k^ 1 ^ 3 with 
wavevector k being in the inertial range, we argue that even for a finite /3, i.e., a finite mean 
magnetic field the energy spectra will scale as k~ 5 ^ 3 in the inertial range, a result supported 
by many observational evidences Q|. The self-consistent amplitude equations will then be 
anisotropic reflecting the presence of underdamped Alfven waves in the inertial range. The 
full correlation matrix will be anisotropic in the hydrodynamic limit; its eigenvalues have 
different amplitudes, but all of them scale the same way. Our confidence on our result that 
the scaling of the correlation functions along directions parallel and perpendicular to the 
direction of the mean magnetic field is same is derived from the fact that our exponent 
values z = 2/3, x — 1/3 satisfy the one-loop integral equations regardless of the strength of 
the mean magnetic field and are consistent with the Ward identity discussed above. 



VI. KOLMOGOROV'S CONSTANTS 

According to the Kolmogorov's hypothesis for fluid turbulence in the inertial range 
energy spectrum E(k) = K Q e 2 ^ 3 k~ 5 ^ 3 , where K Q , a universal constant, is the Kolmogorov's 
constant and e is the energy dissipation rate per unit mass. Various calculations, based 
on different techniques by different groups show that K ~ 1.5 in three 

dimensions. Having noted that the energy spectra, even in the presence of a mean magnetic 
field scale as k~ 5 ^ 3 extensions of Kolmogorov's hypothesis for 3g?MHD allows one to define 
Kolmogorov's constants for the Elsasser fields: E±(k) = K^e± 3 k~ 5 / 3 . Since z ± = u±b, we 
have (z~[(k, uj)z^{— k, —ui)) = (^(k, uj)zJ{— k, —uj)) and e + = e_ = €mhd in absence of any 
crosscorrelations between the velocity and the magnetic fields, we have K„ = K~ = Kmhd- 
The noise strength D and the rate of energy dissipation per unit mass is connected by the 
Novikov's theorem Q, 0:e = 2£>(|^. Noting that the energy spectra E±(k) of z ± in the 
inertial range is given by 

E±(h) = Dk- 3 / V+ k 2 / 3 , (17) 



15 



where D is the effective or renormalised noise strength, we identify 

K MHD = 1.6 [1 + 0.7 (3r 2 - 6r - /3 2 29/105)] 2/3 . (18) 

The notable feature of the expression (|18p is that the constant K depends on the 
dimensionless parameter (3 which we introduced before. For T = 0, i.e., for no magnetic 
fields, we find Kmh d = K Q = 1.6 for pure fluid turbulence which is well-within the accepted 
range of values |3l( . Before closing this Sec. we would like to point out that the presence of 
multiscaling raises questions about Kolmogorov's constant Kmhd being universal: A small 
but finite intermittency correction (i.e., multiscaling) over the simple K41 scaling implies 
the presence of an arbitrary scale which may spoil the universality of Kmhd- We however 
refrain ourselves from getting into this question and adopt a point of view that regardless 
of whether or not Kmhd remains universal due to multiscaling, the numerical value of this 
constant is likely to get affected by the presence of Alfven waves in the system which is 
reflected by the expression (JTSJ). 



VII. POSSIBILITIES OF VARIABLE MULTIFRACTALITY 

Experiments and numerical sumulations find nonlinear multiscaling corrections 

to the K41 prediction of Cp = p/3 for the structure functions in the inertial range. Until the 



account for multiscaling 
assumed a log-normal 



date, no controlled perturbative calculation for C is available, 
in fluid turbulence, however, Obukhov |3J| and Kolmogorov 
distribution for dissipation e to arrive at 

f r \ (*/2)p(p-3) 

S;(r) = (\Av\ p ) = C p I p / 3 r p / 3 (-) , (19) 

where e is the mean value of e and (a bound for e in fluid turbulence has been discussed in 
Ref-laa): 

(e(x + r)e(x)) oc ((At>) 5 /r 2 ) ~ (L/rf 5 . (20) 

For small 5, 5 ~ 95. A standard calculation on the randomly stirred model yields intermit- 
tency exponent S = 0.2 |22j where 5 = 95, whereas the best possible estimate from exper- 
iments is 0.23 j^]. This model, despite having well-known limitations and difficulties 
serves as a qualitative illustration of multiscaling. As in fluid turbulence, in MHD the dissi- 
pations e± of the Elsasser variables fluctuate in space and time, and as a result one may 
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define two intermittency exponents 5± for them. In the present problem S + = <5_ = Smhd- 
Below we calculate the exponent Smhd in a one-loop expression which will give us an esti- 
mate of the Alfven wave speed-dependent deviation of the scaling of the structure functions 
from their simple-scaling values as predicted by the K41 theory. We closely follow Ref.jj^]. 
We work with the self-consistent forms for the self-energies and correlation functions given 
above along with the consistency relations for the amplitude-ratios Tand (3. Following 



Ref. 32], we find the dissipation correlation functions in 3d to be 

(e(x + r)e(x)) ~ 12Ae 2 MHD a 2 K 2 MHD In (21) 

1/3 

with a being defined by the relation u + = ae MHD . From the self-consistent amplitude- 
relations (|13|) we find 

a = 0.4 [1 + 0.7 (3r 2 - 6r - /3 2 29/105)] 1/3 [l - 0.5r - 4/3 2 (0.7 - 0.03r)] 1/3 . (22) 

Thus, a which is a universal coefficient in ordinary fluid turbulence, varies with the parameter 
P, or with the Alfven wave speed in MHD. Substituting the values of K Q and a we find 

S + = 5_ = Smhd = 0.2 [l + 0.7 (3r 2 - 6r - /3 2 29/105)] 4/3 . (23) 

Thus we find that a decrease in the value of S with an increase in (3, i.e., with increasing mean 
magnetic field. Despite the limited applicabilities of log-normal models in characterising 
multiscaling, we can conclude, from our expression ()23|) for the intermittency exponent in 
MHD in the presence of a mean magnetic field, that the intermittency corrections to the 
simple K41 scaling is likely to get affected (reduced in our calculations) in the presence of 
Alfven waves. Further calculations and/or numerical simulations are needed to find the exact 
extent of the dependence of multiscaling on the mean magnetic field and the possibilities of 
anisotropic multiscaling for the structure functions parallel and perpendicular to the mean 
magnetic field as demonstrated recently in Ref-jlT]]. It should be noted that our conclusions 
on the multiscaling properties of the structure functions depend on a log-normal model for 
3<iMHD. Such a description, unfortunately, is unable to distinguish between the structure 
functions parallel and perpendicular to the direction of the mean magnetic field in the 
system. Moreover, the intermittency exponents above [expressions 1)23)1 ] are evaluated in 
the lowest order in mean magnetic field. Therefore, even though from our results we are 
not able to make firm comments on the possibility of parallel and perpendicular structure 
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functions exhibiting different multiscaling, the real importance of our results lie in their 
elucidation of the multiscaling properties depending on the mean magnetic field. 



VIII. CONCLUSION 



In this paper we have considered the effects of the Alfven waves on the statistical proper- 
ties of the correlation functions of the velocity and the magnetic fields or the Elsasser fields. 
We considered the case when the velocity fields are incompressible. In a one-loop approxi- 
mation we find that the effective or the renormalised Alfven wave speed scales as fc _1//3 where 
the wavevector k is in the inertial range. This immediately yields that the energy spectra, 
even in the presence of a mean magnetic field, scale as k~ 5 ^ 3 in the inertial range. We iden- 
tify a dimensionless parameter (3 which is the ratio of the effective Alfven wave speed and 
the renormalised viscosity. We obtain self-consistent relations between the amplitude ratio 
T of the correlation functions and (3. These relations allow us to calculate the dimensionless 
Kolmogorov's constant and we show that it depends explicitly on (3 or on the mean magnetic 
field. Finally, we calculate the intermittency exponent which in a log-normal model gives a 
qualitative account of the multiscaling in terms of the deviation from the K41 scaling for the 
structure functions. We would like to emphasize that although the one-dimensional Burgers- 
like model of MHD of Ref.[2] and the 3<iMHD Eqs. yield energy spectra independent of the 
Alfven waves, the long wavelength physical pictures are different. In the former case, due to 
the nonrenormalization of the bare Alfven wave speed, Alfven waves are overdamped in the 
hydrodynamic limit; the dominant process in that limit is viscous dissipation. In contrast, 
in 3c?MHD, the Alfven wave speed picks up a singular correction in the hydrodynamic limit 
leading to the K41 scaling of the energy spectra and the presence of underdamped Alfven 
waves in the hydrodynamic limit. Despite similar mathematical structures these crucial 
differences arise principally because the one-dimensional model decouples completely when 
written in terms of the Elsasser variables, allowing to comove with the waves of each of them 
separately. As a result correlators of each of them are independent of the Alfven waves and 
the bare Alfven wave speed remains unrenormalised leading to overdamped Alfven waves 
in the hydrodynamic limit. However, the 3dMHD equations do not decouple when written 
in terms of the Elsasser variables and hence oppositely propagating waves cannot be made 
to vanish by comoving. Despite the limitations of the one-loop methods \j\ and the small 
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(5 approximation to facilitate easier analytical manipulations for the amplitude relations, 
we obtain results which are significantly new and open the intriguing possibilities of MHD 
multiscaling universality classes being parametrised by the mean magnetic field. Some of 
the quantitative details will change if one retains terms which are higher order in (3\ however, 
we believe that the qualitative picture will essentially remain the same. As MHD turbulence 
forms a natural example of a driven nonequilibrium system with Alfven waves, our results 
are the first of its kind for a natural system. In our log-normal model approach, we did 
not distinguish between the longitudinal and the transverse structure functions our results 
cannot be compared directly with those of Ref.Q], where the longitudinal and the trans- 
verse structure functions are found to scale differently with multiscaling exponents which 
depend on the magnitude of the mean magnetic field in the presence of a finite magnetic 
helicity. In contrast our results apply to the multiscaling of the usual structure functions 
in the absence of any magnetic helicity which are combinations of the transverse and the 
longitudinal structure functions as considered in Ref. 12] which would also then exhibit a 
mean magnetic field dependent scaling. In particular we find that the usual structure func- 
tions multiscale less in presence of an increasingly strong mean magnetic field which is in 
qualitative agreement with those of Ref. • Although here we have restricted ourselves to 
the study of Alfven waves as relevant perturbations on the amplitude and the scaling of 
the isotropic correlation functions in the limit of small (3, our results indicate the possibil- 
ity of the multiscaling exponents varying with the amplitude of the mean magnetic field. 
Our analyses, when extended for finite values of f3 and for full anisotropic structure of the 
correlation functions, are likely to provide further understanding of and resolve some of the 
discrepancies between the various phenomenological scenarios available for situations with 



large mean magnetic fields 
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, 3]. We leave this task for the future. 
From a broader point of view our results demonstrate the critical role of wave-like ex- 
citations in determining the statistical properties of 3g?MHD. The presence of propagating 
waves is not confined to MHD only; they are a generic feature in many other naturally 
occurring soft-matter systems where such waves can be present in a viscous environment, 



e.g., active polar gels in cytoskeletal dynamics |39j and in the dynamics of self-propelled 



particles 



401 ] etc. We believe our results will lead to similar theoretical and experimental 



studies in relevant nonequilibrium soft-matter systems. 
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